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Abstract Vacuum fluctuations provide a fundamental source of dissipation for sys- 
tems coupled to quantum fields by radiation pressure. In the dynamical Casimir 
effect, accelerating neutral bodies in free space give rise to the emission of real pho- 
tons while experiencing a damping force which plays the role of a radiation reaction 
force. Analog models where non-stationary conditions for the electromagnetic field 
simulate the presence of moving plates are currently under experimental investiga- 
tion. A dissipative force might also appear in the case of uniform relative motion 
between two bodies, thus leading to a new kind of friction mechanism without me- 
chanical contact. In this paper, we review recent advances on the dynamical Casimir 
and non-contact friction effects, highlighting their common physical origin. 



1 Introduction 

The Casimir force discussed in this volume represents the average radiation pressure 
force upon one of the interacting bodies. When the zero-temperature limit is consid- 
ered, the average is taken over the vacuum field state. Although the average electric 
and magnetic fields vanish, the Casimir force is finite because radiation pressure is 
quadratic in the field strength operators. In this sense, the Casimir force derives from 
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the fluctuating fields associated with the field zero-point energy (or more precisely 
from their modification by the interacting bodies). 

As any quantum observable, the radiation pressure itself fluctuates fV,^. For a 
single body at rest in empty space, the average vacuum radiation pressure vanishes 
(for the ground-state field cannot be an energy source), and all that is left is a fluc- 
tuating force driving a quantum Brownian motion |3|. The resulting dynamics is 
characterized by diffusion in phase space, thus leading to decoherence of the body 
center-of-mass |4|. 

Besides diffusion, the radiation pressure coupling also leads to dissipation, with 
the corresponding coefficients connected by the fluctuation-dissipation theorem H. 
As in the classical Brownian motion, the fluctuating force on the body at rest is 
closely related to a dissipative force exerted when the body is set in motion. Since 
the vacuum state is Lorentz invariant, the Casimir dissipative force vanishes in the 
case of uniform motion of a single body in empty space, as expected from the prin- 
ciple of relativity. For a non-relativistic "mirror" in one spatial dimension (ID), 
the Casimir dissipative force is proportional to the second-order derivative of the 
velocity |6|, like the radiation reaction force in classical electrodynamics. Casimir 
dissipation is in fact connected to the emission of photon pairs by the accelerated 
(electrically neutral) mirror, an effect known as the dynamical Casimk effect (DCE). 
The power dissipated in the motion of the mirror is indeed equal to the total radiated 
power in DCE as expected from energy conservation. 

The creation of photons in a ID cavity with one moving mirror was first analyzed 
by Moore |7|, and explicit results were later derived in Ref. |8|. Relativistic results 
for the dissipative Casimir force upon a single mirror in ID and the connection with 
DCE were derived in a seminal paper by Fulling and Davies ||9|- At this earlier 
stage, the main motivation was the analogy with the Hawking radiation associated 
with black-hole evaporation I lOj . 

The interplay between Casimir dissipation and fluctuations was investigated only 
much later lISl fTTIITSI . in connection with a major issue in quantum optics: the fun- 
damental quantum limits of position measurement (this was motivated by the quest 
for interferometric detection of gravitational waves) |13|. Linear response theory 
lfT4l provides a valuable tool for computing the Casimir dissipative force on a mov- 
ing body from the fluctuations of the force on the body at rest, which is in general 
much simpler to calculate. This method was employed to compute the dissipative 
force on a moving, perfectly reflecting sphere |15| and on a plane surface experi- 
encing a time-dependent perturbation |16| (the latter was also computed by taking 
the full time-dependent boundary conditions (b.c.) ifTTlfTSl ). 

In all these configurations, Casimir dissipation turns out to be very small when 
realistic physical parameters are taken into account. The predicted orders of mag- 
nitude are more promising when considering a closely related effect: quantum non- 
contact friction in the shear relative motion between two parallel surfaces |fT9ll2()l . 
In contrast with the radiation reaction dissipative effect discussed so far, quantum 
friction takes place even for a uniform relative motion. On the other hand, quantum 
friction requires the material media to have finite response times (dispersion). From 
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Kramers-Kronig relations ||2T|. the material media must also be dissipative, and the 
resulting friction depends on the imaginary part of the dielectric constant e . 

Whereas the direct measurement of the Casimir radiation reaction dissipative 
force seems to be beyond hope, the corresponding photon emission effect might be 
within reach in the near future. The properties of the radiated photons have been 
analyzed in great detail over recent years. For a single moving plane mirror, the 
frequency spectrum was computed in ID |22J as well as in 3D [23| in the non- 
relativistic approximation. A variety of 3D geometries were considered, including 
deforming mirrors |24|, parallel plates ||25I , cylindrical waveguides f26l, and spher- 
ical cavities containing either scalar ||27]| or electromagnetic fields |28 , 29|. 

Closed rectangular |[30ll3n[32l or cylindrical |33| microwave cavities with one 
moving wall are by far the best candidates for a possible experimental implementa- 
tion, with the mechanical oscillation frequency Q. tuned into parametric resonance 
with cavity field modes. Because of the parametric amplification effect, it is neces- 
sary to go beyond the perturbative approximation in order to compute the intracavity 
photon number even in the non-relativistic regime ll34l [35l [36l . 

As the microwave field builds up inside the cavity, cavity losses (due to trans- 
mission, dissipation or diffraction at the rough cavity walls) become increasingly 
important. Finite transmission at the mirrors of a ID cavity was taken into account 
within the scattering approach developed in Refs. Il22ll37l[38l . Master equations for 
the reduced density operator of the cavity field in lossy 3D cavities were derived in 
Refs. Il39ll40l . Predictions for the total photon number produced at very long times 
obtained from the different models are in disagreement, so that a reliable estimation 
of the DCE magnitude under realistic conditions is still an open theoretical problem. 

It is nevertheless clear that measuring DCE photons is a highly non-trivial chal- 
lenge (see for instance the proposal [41 1 based on superradiance amplification). For 
this reason, in recent years the focus has been re-oriented towards analog models 
of DCE. Although dynamical Casimir photons are in principle emitted even in the 
case of a global 'center-of-mass' oscillation of a cavity, the orders of magnitude 
are clearly more favorable when some cavity length is modulated. In this case, one 
might modulate the optical cavity length by changing the intracavity refractive in- 
dex (or more generally material optical constants) instead of changing the physical 
cavity length. For instance, the conductivity of a semiconductor slab can be rapidly 
changed with the help of a short optical pulse, simulating the motion of a metal- 
lic mirror ll42l |43| and thereby producing photons exactly as in the DCE |44 451. 
An experiment along these lines is currently under way |46| (see Ref. |47| for an 
update). Alternatively, one might select a setup for operation of an optical paramet- 
ric oscillator such that it becomes formally equivalent to DCE with a modulation 
frequency in the optical domain |38 1. 

Alongside the examples in quantum and nonlinear optics, one can also devise 
additional analogues of DCE in the field of circuit quantum electrodynamics BSl . 
For instance, a co-planar waveguide with a superconducting quantum interference 
device (SQUID) at its end is formally equivalent to a ID model for a single mirror 
ll49l l50l . When a time-dependent magnetic flux is applied to the SQUID, it sim- 
ulates the motion of the mirror. More generally, Bose-Einstein condensates also 



4 



Diego A. R. Dalvit, Paulo A. Maia Neto, and Francisco Diego Mazzitelli 



provide interesting analogues for DCE [51] and Casimir-like dissipation ll52l . with 
electromagnetic vacuum fluctuations replaced by zero-point fluctuations of the con- 
densate. 

Reviews on fluctuations and Casimir dissipation on one hand and on DCE on the 
other hand can be found in Refs. Il53ll54ll and Refs. Il55ll56ll57l . respectively. This 
review paper is organized as follows. In Sec. 2, we discuss the interplay between 
fluctuations, dissipation and the photon creation effect for a single mirror in free 
space. Sec. 3 presents a short introduction to non-contact quantum friction. In Sec. 
4, photon creation in resonant cavities with either moving walls or time-dependent 
material properties is presented in detail. Sec. 5 briefly discusses experimental pro- 
posals, and Sec. 6 contains some final remarks. 



2 Dissipative effects of the quantum vacuum 
2.1 ID models 

We start with the simplest theoretical model: a non-relativistic point-like 'mirror' 
coupled to a massless scalar field (j){x,t) in ID. We take the Dirichlet b.c. at the 
instantaneous mirror position q{t) : 

0(^(0,0 = 0. (1) 

In the non-relativistic approximation, we expect the vacuum radiation pressure force 
/(f) to be proportional to some derivative of the mirror's velocity. As a quantum 
effect, the force must also be proportional to h, and then dimensional analysis yields 

/(O^^, (ID) (2) 

where q^"^ (f ) = d"q{t) /dt" . Note that ^ is consistent with the Lorentz invariance of 
the vacuum field state, which excludes friction-like forces proportional to (f ) for 
a single moving mirror (but not in the case of relative motion between two mirrors 
discussed in the next two sections). 

In order to compute the dimensionless prefactor in Q, we solve the b.c. ([T]| to 
first order in q{t) as in Ref. Q, with the mirror's motion treated as a small perturba- 
tion. However, instead of analyzing in the time domain, we switch to the frequency 
domain, which allows us to understand more clearly the region of validity of the 
theoretical model leading to (j2|i. We write the Fourier transform of the field as a 
perturbative expansion: 



(3) 
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where the unperturbed field 0(){x,co) corresponds to a static mirror at x = : 
0o{O,(o) = 0. The boundary condition for 8't>{x,(o) is derived from ([T|i by tak- 
ing a Taylor expansion around x = to first order in Q{£2) (Fourier transform of 

80(0, cOo) = - /~ ^ e(«o - coi) d,<Po{0, CO,)- (4) 

Eq. Q already contains the frequency modulation effect at the origin of Casimir 
dissipation: the motion of the mirror (frequency Q) generates an output amplitude 
at the sideband frequency cOq = COi+ Q proportional to Q{£2) from a given input 
field frequency (Oi . 

In order to find the Casimir dissipative force, we take the Fourier transform of the 
appropriate component Tn — ^[^2(^1^)^ + (<5x0 )^] of the energy-momentum tensor 
and then replace the total field 'P{x, co) containing the solution 8<P{x, co) of eq. (j4]). 
After averaging over the vacuum state, the resulting force is written as 



F{Q)^XmQ{^), (5) 



X{Q)^2i^l -^iQ + (o,)\co,\. (6) 



After regularization, it is simple to show that the contribution jz2 d(0\{...) can- 
cels the contribution f^dcOi{...) in Thus, only field frequencies in the interval 
—Q < (Oi<Q contribute to the dynamical radiation pressure force, yielding 

^^'^) = W2(X2), (7) 

in agreement with ^ with a positive prefactor (g^) as expected for a dissipative 
force. This result was first obtained in Ref. ||6l within the perturbative approach 
and coincides with the non-relativistic limit of the exact result derived much earlier 
in Ref. |9 |. It may also be obtained as a limiting case of the result for a partially 
transmitting mirror, which was derived either from the perturbative approach out- 
lined here [12| or by developing the appropriate Sch winger- Keldysh effective action 
within a functional approach to the dissipative Casimir effect 158 J . When consider- 
ing a moving dielectric half-space, one obtains as expected one-half of the r.h.s. 
of (j?]) in the limit of an infinite refractive index f59'|. The final result for F{Q) is 
exactly the same as in (j?]) when we replace the Dirichlet b.c. ([T]) by the Neumann 
b.c. at the instantaneously co-moving Lorentz frame (primed quantities refer to the 
co-moving frame) <3y0 |v'=^(f') = [60|. Dirichlet and Neumann b.c. also yield the 
same force in the more general relativistic regime 1,61.1 . On the other hand, for the 
Robin b.c. ^ 

dx'(^\x'=q{t') ^ -p(l'lx'=q{t') (8) 

(/3 is a constant parameter), the force susceptibility x(i2) displays a non-monotonic 
dependence on Q, and is nearly suppressed at X2 ~ 2.5c/ P 
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We could have anticipated that high-frequency modes with Oi S> i2 would not 
contribute because they "see" the mirror's motion at frequency 12 as a quasi- 
static perturbation, and indeed the dissipative Casimir force originates from low- 
frequency modes for which the motion is non-adiabatic. But there is a more illu- 
minating interpretation that explains why the contribution comes precisely from the 
frequency interval [—12,0]. In Fig.[r| we show that this interval corresponds to the 
field modes leading to frequency sidebands cOo ~ COi+ £2 (see eq. Q) across the 
border between positive and negative frequencies (for Q <0 the corresponding in- 
terval is [0, —Q] and the analysis is essentially the same). In other words, for these 
specific modes the motional frequency modulation leads to mixing between posi- 
tive and negative frequencies. Bearing in mind the correspondence between positive 
(negative) frequencies and annihilation (creation) operators, this mixing translates 
into a Bogoliubov transformation coupling output annihilation operators to input 
creation ones, and viceversa [63] (examples will be presented in Section 4). The 
important conclusion is that sideband generation for these modes corresponds to 
photon creation (and also annihilation in the case of a general initial field state), 
whereas outside the interval [—12,0], where no mixing occurs, the sideband effect 
corresponds to photon inelastic scattering with neither creation nor annihilation. 

From this discussion, we can also surmise the important property that the dy- 
namical Casimir photons have frequencies bounded by the mechanical frequency 
Q, as long as the perturbative non-relativistic approximation holds |22|. More- 
over, the Casimir photons are emitted in pairs, with photon frequencies satisfying 
0)1 + 0)2 = Q. Hence the frequency spectrum is symmetrical around Q /2, where it 
has a maximum for the Dirichlet case 122] but not generally for the Robin b.c. Il64ll . 

In short, the derivation of (j2]) from (j4]) highlights the direct connection between 
the dissipative dynamical force and the dynamical Casimir photon emission effect. 
The dissipative force thus plays the role of a Casimir radiation reaction force, damp- 
ing the mirror's mechanical energy as Casimir photons are emitted. In fact, the ex- 
pression given by (|2]) has the same form as the radiation reaction force in classical 



I 
I 

CO, |cO|+Q 
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-Q frequency 



Fig. 1 For a given mechanical modulation frequency £2 (assumed to be positive in this diagram), 
the field modes contributing to the dissipative dynamical Casimir effect lie in the interval [— 12 , 0] , 
which corresponds to negative input field frequencies yielding positive sideband frequencies. 
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electrodynamics, apart from a dimensionless pre-factor inversely proportional to the 
fine-structure constant /{he). 

The discussion on the field modes actually contributing in (j4]) also allows us to 
address the domain of validity of the various assumptions employed in the derivation 
presented above. We have assumed in ([T]) that the field vanishes at the instantaneous 
mirror's position, no matter how fast the mirror and field oscillate. However, the 
electric currents and charge density inside a real metallic mirror respond to field 
and position changes over a finite time scale, so that we expect our oversimplified 
model to be physically meaningful only at low frequencies (more general results 
and discussions are presented in Refs. 1 12, 65 1). Since the relevant field frequencies 
are bounded by the mechanical frequencies {\(0\\ < \Q |), the model is consistent as 
long as typical mechanical frequencies are much smaller than the frequency scales 
characterizing the metallic response - typically the plasma frequency of metals. 

A second point to be clarified is the connection between the perturbative linear 
approximation employed above and the non-relativistivistic approximation. When 
deriving Q from ([T]i, we have taken the long-wavelength approximation to expand 
the field around x = Q. Let us assume, to simplify the discussion, that the mirror os- 
cillates with frequency Q and amplitude qo. The non-relativistic regime then trans- 
lates into Qqo/c 1 and all relevant field modes correspond to long wavelengths 
2nc/cOi ^ ^0 since they satisfy the inequality \cOi\ < Q. More generally, the long- 
wavelength approximation follows from the non-relativistic condition provided that 
there is an inertial reference frame for which the motion is spatially bounded. 



2.2 Casimir-driven decoherence 

As in classical Brownian motion, the dissipative effect is closely related to the fluc- 
tuations provided by the reservoir, which in our case is the quantum vacuum field. 
This connection provides yet another tool for computing the dissipative Casimir 
force: by using linear response theory, one can derive (j7]l from the correlation func- 
tion of the force on a static plate lfT2l . (This method has been employed for different 
geometries in 3D [ Tn[T5l[T6ll26ll ). More interestingly, if we take the mirror's posi- 
tion as a full dynamical observable rather than a prescribed function of time, vacuum 
radiation pressure plays the role of a Langevin fluctuating force Il3]|58j, leading to 
diffusion in phase space, which adds to the associated average dissipative force (j7]i. 

Let us first analyze the mechanical effect of the dissipative Casimir force in the 
context of classical dynamics. We consider a point-like mirror of mass M in a har- 
monic potential well of frequency Q . Taking the Casimir dissipative force given by 
(jTj) into account, and neglecting for the time being any associated stochastic force, 
we write the mirror's equation of motion as 
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We assume that the oscillator's zero-point energy hQ/2 is much smaller than the 
rest mass energy Mc^, and then find oscillatory solutions of ^ which are damped 
at the rate 

1 hQ 

This result provides a good illustration of how weak the Casimir dissipation is for a 
single mirror in a vacuum. On the other hand, the associated diffusion in phase space 
is more relevant, particularly in the context of the full quantum theory discussed be- 
low, since it provides an efficient decoherence mechanism for non-classical quantum 
states. 

The quantum description of the mirror's dynamics can be developed from the 
Hamiltonian for the radiation pressure coupling with a dispersive semi-transparent 
mirror (transparency frequency (Oc) ll66l. One derives a master equation for the re- 
duced density operator p of the mirror, which can also be cast in the form of a 
Fokker-Planck equation for the Wigner function W{x,p,t) representing the mirror 
quantum state [4J : 

d,W = -(I - AM/M)—d,W +Mnhd„W + 2rdJxW)+Di^-^ - D2^r^- 
M ' ' dx^ dxdp 

(11) 

The time-dependent coefficients AM (mass correction), F (damping coefficient), 
Di and D2 (diffusion coefficients) are written in terms of the correlation function 
of the field linear momentum operator. The perfectly reflecting limit corresponding 
to ^ is obtained when i2 ^ 0)^, in line with our previous discussion since l/(Oc 
represents the characteristic time scale for the material medium response. In this 
limit, for t>l/ cOc, r{t) rapidly approaches the expected constant value as given by 



(10 1, whereas £>i (f ) approaches the asymptotic value 

Di=hr/{MQ) (12) 

for t >2n/Q. This connection between diffusion and damping plays the role of a 
fluctuation-dissipation theorem for the vacuum state (zero temperature). 

Under the time evolution described by the Fokker-Planck equation (111, an ini- 
tially pure state gradually evolves into a statistical mixture. The physical reason be- 
hind this decoherence effect is the buildup of entanglement between the mirror and 
the field due to the radiation pressure coupling and the associated dynamical Casimir 
photon creation |4|. As an example of initial state, we consider the "Schrodinger's 
cat" superposition of coherent states 1 1//) = (|a) + | — a))/\/2 with the amplitude 
a = iPo/^/2MhcOQ along the imaginary axis (zLPq are the average momenta corre- 
sponding to each state's component). The corresponding Wigner function W{x,p) 
contains an oscillating term proportional to cos(2Po-^/^) which represents the coher- 



ence of the superposition. Clearly, the diffusion term proportional to Di in ( 11 1 will 



wash out these oscillations along the position axis, thus transforming the cat state 



into a mixture of the two coherent states. From (111 and ( 12 1, the corresponding 
decoherence time scale is 
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t., = J^-^(^] r-S (13) 



IP^Di \2Po 

where dp — ^ MhQ /2 is the momentum uncertainty of the coherent state (we have 
assumed that Qt^ S> 1). Since 2Po measures the distance between the two compo- 



nents in phase space, ( 13 i shows that decoherence is stronger when the state compo- 
nents are further apart, corresponding to thinner interference fringes in phase space. 
Decoherence from entanglement driven by dynamical Casimir photon creation is 
thus very effective for macroscopic superpositions in spite of the smallness of the 
corresponding damping coefficient F. As the radiation pressure control of microres- 
onators improves all the way to the quantum level ||67]| . Casimir driven decoherence 
might eventually become of experimental relevance. 



2.3 3D models 

The orders of magnitude can be reliably assessed only by considering the real-world 
three-dimensional space. We start with the simplest geometry in 3D: a plane mirror 
parallel to the xy plane, of area A and moving along the z— axis. For an infinite plane, 
we expect the dissipative Casimir force to be proportional to A, so that we have to 
modify (|2| to include a squared length. For a scalar field satisfying a Dirichlet b.c. 
analogous to ([T]l, the derivation is very similar to the one outhned above (S): 

= - ^^^1^^) scalar) (14) 

For the electromagnetic field, the model of perfect reflectivity provides an accu- 
rate description of metallic mirrors at low frequencies. For a mirror moving along 
the z-axis, the electric and magnetic fields in the instantaneously co-moving Lorentz 
frame S' satisfy 

ZXE'|n,iiTor = 0, Z-B'Imin-or = 0. (15) 

It is useful to decompose the fields into transverse electric (TE) and transverse 
magnetic (TM) polarizations (where 'transverse' means perpendicular to the inci- 
dence plane defined by z and the propagation direction). For the TE component, we 
define the vector potential in the usual way: E'™^ = — (9fA(™\ B'^^' = V x A'^^' 
under the Coulomb gauge V • A^^^' = 0. Since A*^^^^ • z = 0, A^™' is invariant under 
the Lorentz boost from the co-moving frame to the laboratory frame. The resulting 
b.c. is then similar to ([T]): 

A^^''\x,y,q{t),t)^0. (16) 

On the other hand. A'™' has a component along the z-axis, so that the Coulomb 
gauge is no longer invariant under the Lorentz boost, resulting in complicated b.c. 
also involving the scalar potential. It is then convenient to define a new vector po- 
tential j^/'^™^ as 
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E(™) = V X ^(™) ^ g(TM) ^ ^^^(TM) 

under the gauge = 0. Like A'^^', is also invariant under Lorentz 

boosts along the z-axis. From ( 15 i, one derives that satisfies a Neumann b.c. 



at the co-moving frame, yielding 

[d,+q{t)dr + ffii)] J^^™\x,y,q{t),t) = 0. (17) 

The condition of perfect reflectivity then results in two independent problems: a 
Dirichlet b.c. for TE modes, and a Neumann b.c. in the instantaneously co-moving 
frame for TM modes. The TM contribution turns out to be 11 times larger than the 
TE one, which coincides with ( [T4| . The resulting dissipative force is then ll68l 

Mfl(^'ff) 

fit) = - J 2 4 ■ (3D, electromagnetic) (18) 

As in the ID case, the dissipative Casimir force plays the role of a radiation reac- 
tion force, associated with the emission of photon pairs with wave-vectors satisfying 
the conditions |ki| + |k2| =i2/candki|| =— k2|| from translational symmetry paral- 
lel to the plane of the mirror The angular distribution of emitted photons displays an 
interesting correlation with polarization: TE photons are preferentially emitted near 
the normal direction, whereas TM ones are preferentially emitted at larger angles, 
near a grazing direction if the frequency is smaller than Q /2 ll23l . 

Results beyond the model of perfect reflection were obtained in Ref |69| for a 
dielectric half-space (see also Ref. |70| for the angular and frequency spectra of 
emitted photons). In this case, there is also photon emission (and the associated 
dissipative radiation reaction force) if the dielectric mirror moves sideways, or if a 
dielectric sphere rotates around a diameter |71 1. We will come back to this type of 
arrangement when discussing non-contact quantum friction in the next section. 

To conclude this section, we compute the total photon production rate for a per- 



fectly reflecting oscillating mirror directly from ( 18 i. By energy conservation, the 
total radiated energy is the negative of the work done on the mirror by the dissipative 
Casimir force: ^ 

E = - r f{t)q'^'Ht)dt. (19) 



We evaluate the integral in ( 19 1 using the result ( [T8] l for an oscillatory motion of 
frequency Q and amplitude go exponentially damped over a time scale T ^ l/Q : 

Since the spectrum is symmetrical with respect to the 
frequency i2/2, we can derive the number of photons from the radiated energy 
using the relation E = NhQ /2. The total photon production rate is then given by 



Fluctuations, dissipation and the dynamical Casimir effect 



11 



with Vmax = •^2<7o and Ao = 2nc/ Q. representing the typical scale of the relevant 
wavelengths. With Vmax/c ^ 10"^, Q./{2n) ^ 10 GHz and A ^ A.^ ^ lOcm^, we 
find N/T ^ 10^^ photons/sec or approximately one photon pair every two days! 

The dynamical Casimir effect is clearly very small for a single oscillating mirror. 
Adding a second parallel plane mirror, the photon production rate is enhanced by a 
factor {QL/c)^^ ^10^ for separation distances L in the sub-micrometer range Il25l . 
But at such short distances, finite conductivity of the metallic plates, not considered 
so far, is likely to reduce the photon production rate. In this type of arrangement, a 
much larger effect is obtained by considering the shear motion of one plate relative 
to the other (instead of a relative motion along the normal direction). As discussed 
in the next section, because of finite conductivity a large friction force is predicted 
at short distances, which results from the creation of pairs of excitations inside the 
metallic medium 1 20 , 72 1. As for the emission of photon pairs, the orders of magni- 
tude are more promising when considering a closed cavity with moving walls, to be 
discussed in Section 4. 



3 Quantum Friction 

There is an intimate connection between the dynamical Casimir effect and the possi- 
bility that electrically neutral bodies in relative motion may experience non-contact 
friction due to quantum vacuum fluctuations, the so-called "quantum friction". As 
we have discussed in the previous section, dielectric bodies in accelerated motion ra- 
diate Casimir photons. Shear motion of two bodies, even at constant relative speed, 
can also radiate energy. Just as in the case of a single accelerated mirror in a vac- 
uum, shear motion cannot be removed by a change of reference frame. A frictional 
force between two perfectly smooth parallel planes shearing against each other with 
a relative velocity v results from the exchange of photons between the two surfaces. 
These photons carry the information of the motion of one surface to the other one, 
and as a result linear momentum is exchanged between the plates, leading to friction. 

In order to illustrate the physics of quantum friction we will follow here an ap- 
proach due to Pendry [20| who considered the simplest case of zero temperature 
and the non-retarded (van der Waals) limit. The nice feature of this approach is that 
it manifestly connects to the intuitive picture of motion-induced (virtual) photons as 
mediators of momentum exchange between the shear surfaces. A dielectric surface, 
although electrically neutral, experiences quantum charge fluctuations, and these 
have corresponding images on the opposing dielectric surface. Since the surfaces 
are in relative parallel motion, the image lags behind the fluctuating charge distribu- 
tion that created it, and this results in a frictional van der Waals force. Note that for 
ideal perfect metals, the image charges arrange themselves instantaneously (do not 
lag behind), and therefore no quantum friction is expected in this case. 

We model each of the dielectric surfaces as a continuum of oscillators with 
Hamiltonian 
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Ha = 2^^<*^a;k7(fla.k/a;k> + 1/2), (21) 

where aj^ ty and aa;k/ are creation and annihilation bosonic operators associated 
with the upper (a = m) or lower (a = /) plate. Each mode on each surface is defined 
by k, which is a wave-vector parallel to the planar surface, and by which denotes 
degrees of freedom perpendicular to the surface. Following Pendry we restrict our- 
selves to the non-retarded limit (very long wavelengths for the EM modes). In this 
limit the EM field is mainly electrostatic, only the static TM polarization matters for 
a dielectric surface (since the static TE field is essentially a magnetic field that does 
not interact with the non-magnetic surface), and the intensity decays exponentially 
from the surfaces (evanescent fields). The coupling between the oscillator modes 
belonging to different surfaces is mediated by the EM field, and it is assumed to be 
a position-position interaction of the form 

i/i„t (f ) = E <^ki/ {d)x,,^j ® i,._k/ e-'^'"" ■ (22) 
k;/ 

In the non-relativistic limit the effect of the surfaces shearing with speed v along the 
X direction is contained in the last exponential factor. This type of Hamiltonian fol- 
lows from the effective electrostatic interaction between the fluctuating charges in 
the dielectrics and its expansion to lowest order in the displacement of each oscilla- 
tor from its equilibrium position (equivalently, it also follows from the non-retarded 
and static limit of the dipole-dipole interactions between fluctuating dipoles in each 
surface). The coupling factors C]^jji{d) can be obtained by analyzing how each os- 
cillator dissipates energy into the vacuum gap. This is done in Ref. Il20l in two ways, 
by invoking a scattering type of approach relating the fields at the interphases with 
reflection amplitudes, and by considering how the fluctuating charge distributions in 
each surface dissipate energy. The result is C]ijf{d) — (j3k)j3ky/2^eo)e^l''l'', where 



e{co)-l 
e(to) + l 



(23) 



Note the exponential decay due to the evanescent nature of the EM field. In this 
equation dN/dO) is the density of oscillator modes at frequency o and e(ft)) is the 
complex dielectric permittivity of the plates (assumed to be identical). Although 
the interaction Hamiltonian does not depend explicitly on the quantized EM field 
(because this derivation is semiclassical), one can infer the quantum processes of 
creation and absorption that take place by expanding the product Xy tj (8ii/;_k/- 

■^«;kj = - ««,k>] ® [4-k/ ^ «/.-k/] • (24) 

Imagine the system of the two dielectrics is initially in the ground state at zero 
temperature, \\^{t = 0)) = \\irg)u®\Wg)i, where |v/g)a Ilk; I ^s-.k;) a is the prod- 
uct of the harmonic oscillators' ground states for surface a. Eq.([24| implies that 
two motion-induced virtual photons created from an EM vacuum produce one 
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excitation in each surface, i.e., there is a non-zero probabihty of transition to 
states \ l;kj)uiS) |1;— k/)/. The transition probability can be computed using time- 
dependent perturbation theory for the perturbation Hi^tit)- To first order, the transi- 
tion probability from the ground state into each of these two-excitation states is 



4 sin^ [{co„,kj + coi^^kf - kxv)t/2] 



(25) 



In the limit of large times (f oo) we use that sin^{Qt /2)/{Q jT)^ w ntb{Q.) (here 
5(i2) is Dirac's delta function), and therefore the transition probability grows lin- 
early in time. We can find the frictional force equating the frictional work F^v with 
the rate of change in time of the energy of the excitations, namely 



dU 



/i(co„,kj + «/-k/)- 



dt 



(26) 



and the r.h.s. is time-independent since the transition probabilities grow linearly 
in time. Using the expression (23i for j3^y and the transition probabilities (25i at 
large times, and writing the sums over the dielectric degree of freedom j as 27, 

d(i)dN /d(£) (and similarly for /), one finally obtains the following expression for 
the frictional force 



d^k 



.-21kW 



l-kxV 

Jo ' 



'e{co)-l' 


Im 


_e(a)) + l_ 





£{kxV-CO) - 1 
£{kjV-CO) + 1 



(27) 



In the literature there are other more rigorous approaches to calculate quantum 
friction that go beyond the non-retarded quasi-static limit considered above, and 
that can take into account effects of relativistic motion as well as finite temper- 
ature. One of these approaches [73 1 follows the spirit of the Lifshitz-Rytov the- 
ory ll74l . considering the fluctuating electromagnetic (EM) field as a classical field 
whose stochastic fluctuations satisfy the fluctuation-dissipation relation that relates 
the field fluctuations with the absorptive part of the dielectric response of the plates. 
The EM field is a solution to Maxwell's equations with classical fluctuating current 
densities on the plates as source fields, and it satisfies the usual EM b.c. imposed 
on the comoving reference frames on each plate. The relation between the EM 
fields in different frames is obtained via Lorentz transformations. An alternative 
full quantum-mechanical approach considers the quantum EM field in interaction 
with (quantized) noise polarizations and noise currents within the plates |75|. As 
before, the fields in each reference frame are related by Lorentz transformations. In 
this approach the quantum expectation value of the noise currents is given by the 
(quantum) fluctuation-dissipation relation. 

Quantum friction can also happen for neutral atoms moving close to surfaces. 
The theoretical methods to compute the frictional force in these cases are similar to 
the surface-surface quantum friction, and we refer the reader to some of the relevant 
works Il76l 17711781. 
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4 Resonant photon creation in time dependent cavities 

As mentioned in the Introduction, photon creation can be enhanced in closed cavi- 
ties: if the external time dependence involves a frequency that is twice the frequency 
of a mode of the electromagnetic field in the (unperturbed) cavity, parametric am- 
plification produces a large number of photons. As we will see, under certain cir- 
cumstances (ideal three dimensional cavities with non equidistant frequencies in the 
spectrum) the number of photons in the resonant mode may grow exponentially. 
Parametric amplification can take place by changing the length of the cavity with 
a moving surface, but also by changing its effective length through time dependent 
electromagnetic properties of the cavity. 

In order to simplify the notation, in this section we will use the natural units 



4.1 Dynamical Casimir effect in ID cavities 

As in Section 2, we start with a massless real scalar field in a ID cavity with one 
mirror fixed at x = and the other performing an oscillatory motion 



where Q. is the external frequency and e ^ 1 . As we will be mainly concerned with 
situations where QLq ~ 0{\), the maximum velocity of the mirror will be of order 
e, and therefore small values of e correspond to a non-relativistic motion of the 
mirror. We shall assume that the oscillations begin at f = 0, end ait ~T, and that 
L{t) — Lq for f < and t > T. The scalar field (j){x,t) satisfies the wave equation 
□0 = and Dirichlet b.c. 0(x = 0,f) — ^{L{t),t) — 0. When the mirror is at rest, 
the eigenfrequencies are multiples of the fundamental frequency tt/Lq. Therefore, in 

order to analyze resonant situations we will assume that Q. — qn/LQ, q— 1,2,3, 

Inside the cavity we can write 



where the mode functions \j/ii{x,t) are positive frequency modes for f < 0, and ai^ and 
are time-independent bosonic annihilation and creation operators, respectively. 
The field equation is automatically verified by writing the modes in terms of Moore's 
function 7J(f) Q as 



L{t) = Lq[1 + esm{nt)] , 



(28) 



(29) 



k=l 




(30) 
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The Dirichlet boundary condition is satisfied provided that R{t) satisfies Moore's 
equation 

R{t + L{t))~R{t~L{t))^2. (31) 

These simple expressions for the modes of the field are due to conformal invariance, 
a symmetry for massless fields in one spatial dimension. 

The solution to the problem involves finding a solution R{t) in terms of the pre- 
scribed motion L{t). For f < the positive frequency modes are given hy R{t)=t/LQ 
for —Lq <t< Lq, which is indeed a solution to Eq.( 3 1 1 for f < 0. For f > 0, Eq.([3T]l 
can be solved, for example, using a perturbative expansion in e similar to the one 
employed in Section 2 for a single mirror. However, as the external frequency is 
tuned with the unperturbed modes of the cavity, in general there will be resonant 
effects, which produce secular terms proportional to e"'{Qt)" with m < n. Thus 
this expansion is valid only for short times eQt ^ 1. It is possible to obtain a non- 



perturbative solution of Eq.(31 1 using Renormalization Group (RG) techniques ||34|| . 
The RG-improved solution automatically adds the most secular terms, [eQt)", to all 
orders in £, and is valid for longer times e^Q.t <^ 1. The RG-improved solution is 

ElCSl 

t 2 ' 



m = r- 



Kq 



-Imln 



(32) 



25.0 




Fig. 2 R(t) vs. //Lo as given by Eq.|32|. The values of the parameters are q — \ and e = 0.01. 



16 



Diego A. R. Dalvit, Paulo A. Maia Neto, and Francisco Diego Mazzitelli 



where ^ — exp[(— l)'f+'7rg'ef/Lo]- As shown in Fig|2] the function R{t) develops a 
staircase shape for long times ||34l |80l . Within regions of f between odd multiples 
of Lo there appear ^ jumps, located at values of f satisfying cos{q7Zt/Lo) ~ =f1, the 
upper sign corresponding to even values of q and the lower one to odd values of q. 

The vacuum expectation value of the energy density of the field is given by ||9] 
{Too{x,t)) = -/(f+x) -/(f-x), where 



R' 2\R' 2 



(33) 



For q—l ("semi-resonant" case) no exponential amplification of the energy density 
is obtained, whereas for ^ > 2 ("resonant" cases) the energy density grows exponen- 
tially in the form of q traveling wave packets which become narrower and higher as 
time increases (see Fig|4]). Note that, as the energy density involves the derivatives 
of the function R{t), there is one peak for each jump of R{t). 

The number of created particles can be computed from the solution given by 



Eq.( 32 1. Photons are created resonantly in all modes with n = q + 2j, with j a non- 
negative integer This is due to the fact that the spectrum of a one dimensional 
cavity is equidistant: although the external frequency resonates with a particular 
eigenmode of the cavity, intermode coupling produces resonant creation in the other 
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Fig. 3 Energy density profile between plates for fixed time t/Lo = 20.4 for the Cj = 4 case. The 
amplitude coefficient is e = 0.01 



Fluctuations, dissipation and the dynamical Casimir effect 17 

modes. At long times, the number of photons in each mode grows linearly in time, 
while the total number of photons grows quadratically and the total energy inside 
the cavity grows exponentially |30|. These different behaviors are due to the fact 
that the number of excited modes, i.e. the number of modes that reach a growth 
linear in time, increases exponentially. 

The production of massless particles in one dimensional cavities has been ana- 
lyzed numerically in Ref. |81 1. As expected, the numerical evaluations are in perfect 
agreement with the analytic results described above in the case of small amplitudes 
£ < 0.01. 

The force on the moving mirror can be computed as the discontinuity of (Tn) at 
X — L{t). The force produced by the outside field is one -half the expression derived 
in Section 2. This is much smaller than the intracavity contribution and will be 
neglected. Therefore 

(F)«(r„(L(f),f)) = (roo(L(f),f)), (34) 

where the energy momentum tensor is evaluated inside the cavity. This expression 
reproduces the usual attractive result when the miiTor is at rest (f < 0) 

^^^ = -2^1' '''' 

that is, the static Casimir effect in ID. However, at long times it becomes an expo- 
nentially increasing pressure due to the presence of real photons in the cavity ll30l . 

All this treatment can be extended to the case of Neumann b.c. n''<3^0 1 mirror = 0^ 
where is a unit two-vector perpendicular to the trajectory of the mirror. The 
modes of the field can be written in terms of Moore's function R{t) as 



Note the change of sign between Eq.( 30 1 for Dirichlet modes, and Eq.( 36 1 for Neu- 
mann modes. The spectrum of motion-induced photons is the same for both Dirich- 
let and Neumann b.c. |60|, but not for the mixed configuration with one Dirichlet 
mirror and one Neumann mirror 1821 . The one dimensional DCE has also been in- 
vestigated for cavities with Robin b.c. 



4.2 Photon creation in 3D cavities 



The one dimensional DCE with Dirichlet b.c. described in the previous section is 
not only of academic interest: it describes photon creation for the TEM modes in a 
3D cylindrical cavity with a non-simply connected section lf33l . However, in order 
to analyze TE and TM modes in general 3D cavities, a new approach is needed, 
since conformal invariance is no longer useful in 3D. 
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We shall first describe in some detail the simpler case of a scalar field in a rect- 
angular cavity satisfying Dirichlet b.c. |[3T| . and then comment on the extension to 
the case of the electromagnetic field in cylindrical cavities with an arbitrary section. 



4.2.1 Scalar field 



We consider a rectangular cavity formed by perfectly reflecting walls with dimen- 
sions Lx,Ly, and L,. The wall placed at z = L, is at rest for t < and begins to 
move following a given trajectory, L^{t), at f = 0. We assume this trajectory as 
prescribed for the problem (not a dynamical variable) and that it works as a time- 
dependent boundary condition for the field. The field (j) (x, t) satisfies the wave equa- 
tion 0(j) ~ 0, and the b.c. [walls — for all times. The Fourier expansion of the field 
for an arbitrary moment of time, in terms of creation and annihilation operators, can 
be written as 

^ix,t)=Y,a'^u„{x,t)+}i.c., (37) 

n 

where the mode functions Ua{x,t) form a complete orthonormal set of solutions of 
the wave equation with vanishing b.c. 

Wlien f < (static cavity) each field mode is determined by three positive integers 
rixjity and n^, namely 



Mn(x,r < 0) 




(38) 



with (On = n^ (g) n- yj^^j , yj^: 

When f > the boundary condition on the moving wall becomes (j){x,y,z 



L^{t),t) — 0. In order to satisfy it we expand the mode functions in Eq.(37i with 
respect to an instantaneous basis 
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<2r'(o) 



eL"^(o) 



(40) 



In this way we ensure that, as long as L^{t) and L^{t) are continuous at f = 0, each 
field mode and its time derivative are also continuous fiinctions. The expansion in 
Eq.([39]l for the field modes must be a solution of the wave equation. Taking into 
account that the (jOk's form a complete and orthonormal set and that they depend on 
t only through L,(f), we obtain a set of (exact) coupled equations for Qm^ (f) ||3T| : 



Qt^ + <(0 gL"^ = 21 it) ^g^j e]"' + A (f ) 



l(n) 



A'(f)i:§im^ije]"', 



(41) 



where 



(Om{t) = It 




4(0 



(42) 



The coefficients gj^j are defined by 



dz 



d(p„ 



(Pi- 



rn 



The annihilation and creation operators aJJ, and correspond to the particle 
notion in the 'in' region (t < 0). If the wall stops for f > fflnai, we can define a new 
set of operators, a™' and dj^"\ associated with the particle notion in the 'out' region 
{t > ffinai)- These two sets of operators are connected by means of the Bogoliubov 
transformation 

< -E(«;;«nm+4"i3*™). (44) 



The coefficients anm and jSnm can be obtained as follows. When the wall returns to 



its initial position the right hand side in Eq.(41 1 vanishes and the solution is 



6m' (f > ffinal) 



(45) 



with and being some constant coefficients to be determined by the con- 
tinuity conditions at f = fflnai- Inserting Eq.(45i into Eqs.(37i and (39i we obtain 
an expansion of (j) in terms of aJJ, and a^" for t > tf[„^\. Comparing this with the 
equivalent expansion in terms of aj^' and am"' it is easy to see that 



an 



(n) 



(n) 
m • 



(46) 



The amount of photons created in the mode m is the average value of the number 
operator am"'ai^' with respect to the initial vacuum state (defined through am|Oin) — 
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0). With the help of Eq.(|44]l and Eq.(|46]l we find 

i^m) = (0.„ I fll™'«™' I 0,„) ^ Y. 2®m lA'^m^l^- (47) 

n 

In the approach described so far we worked at the level of the dynamical equation 
for the quantum scalar field. Alternatively, one can analyze the problem using the 
effective Hamiltonian method developed in Ref. [85|. The idea is the following. 
Assume that a massless scalar field is confined within a time dependent volume and 
satisfies Dirichlet b.c. At the classical level, the field can be expanded in terms of a 
basis of functions /a(x,r) that fulfill the b.c. at each time, that is 

0(x,f)=^<?„(f)/„(x,f). (48) 

a 

For the rectangular cavities considered in this section these functions can be cho- 
sen to be (f)fa(^,L^{t)). Inserting this expansion into the Klein-Gordon Lagrangian, 
one ends up with a Lagrangian for the generalized coordinates qa{t), which is a 
quadratic function of qa{t) and ^c((f), i.e. it describes a set of coupled harmonic os- 
cillators with time dependent frequencies and couplings. This system can be quan- 
tized following the usual procedure, and the final results for the number of created 
photons are equivalent to those obtained in Eq. ( [47] l. 



4.2.2 Parametric amplification in 3D 



As in Section 4.1, we are interested in resonant situations where the number of pho- 
tons created inside the cavity could be enhanced for some specific external frequen- 



cies. So we study the trajectory given in Eq.(28 1. To first order in £, the equations 



for the modes Eq.(41 1 take the form 



~l7 



sm{Qt)Q^^^ ~ eQhiniQt) Y,gmjQ\ 



(n) 



+ 2ef2cos(i2f)^g™je]"' + 0(e2 
j 



(49) 



It is known that a naive perturbative solution of these equations in powers of 
the displacement e breaks down after a short amount of time, of order {e£2)^^. As 
in the ID case discussed in the previous section, this happens for those particular 
values of the external frequency Q such that there is a resonant coupling with the 
eigenfrequencies of the static cavity. In this situation, to find a solution valid for 
longer times (of order e^^i2^') we proceed as follows. We assume that the solution 



of Eq.(49 1 is of the form 



Gm'(f)=AL"'(fy 



(50) 
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where the functions and B\a are slowly varying. In order to obtain differen- 
tial equations for them, we insert this ansatz into Eq.([49]l and neglect second-order 
derivatives of A^^^ and B^' • After multiplying the equation by e*"*^' we average 
over the fast oscillations. The resulting equations are 



1 dA, 



(n) 



e dt 



L 

j 



L(-®j + y)^(-®m-®j + ^)^ 



Q Q 



(51) 



and 



1 dB, 



L 

j 



2av 



(52) 



where we used the notation 8{co) for the Kronecker 5-function 5a,o- 

The method used to derive these equations is equivalent to the "multiple scale 
analysis" |86| and to the slowly varying envelope approximation ll87l . The equa- 
tions are non-trivial (i.e., lead to resonant behavior) if i2 = IcOm (resonant condi- 
tion). Moreover, there is intermode coupling between modes j and m if any of the 
conditions | tOm ± tOj | = 12 is satisfied. 

We derived the equations for three dimensional cavities. It is easy to obtain the 
corresponding ones for one dimensional cavities. The conditions for resonance and 
intermode coupling are the same. The main difference is that for one dimensional 
cavities the spectrum is equidistant. Therefore an infinite set of modes may be 
coupled. For example, when the external frequency is i2 = 2(0i, the mode m is 
coupled with the modes m±2. This has been extensively studied in the literature 



In what follows we will be concerned with cavities with non-equidistant spec- 
trum. Eqs.([5T| and (|52]) present different kinds of solutions depending both on the 
mirror's frequency and the spectrum of the static cavity. In the simplest 'parametric 
resonance case' the frequency of the mirror is twice the frequency of some unper- 
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turbed mode, say £2 — 2a)ni- In order to find A^'' and B^' from Eq.(|51| and Eq.(|52 1 



we have to analyze whether the coupHng conditions | (Om ± Wj | = ^2 can be satisfied 
or not. If we set Q. = IcOm, the resonant mode m will be coupled to some other mode 
j only if Oj — Oni = £2 ~ 2(£)m- Clearly, the latter relation will be satisfied depending 
on the spectrum of the particular cavity under consideration. 

Let us assume that this condition is not fullfilled. In this case, the equations for 
can be easily solved and give 



A^^ and B 



i^m) = sinh^ 



~l7 



(53) 



In this uncoupled resonance case the average number of created photons in the mode 
m increases exponentially in time. Another way of looking at this particular situation 
is to note that, neglecting the intermode couplings, the amplitude of the resonant 
mode satisfies the equation of an harmonic oscillator with time dependent frequency. 
For the particular trajectory given in Eq.([28]l, the dynamics of the mode is governed 
by Eq. (|49]l with g^j = , that is 



(n) 



coi-2£ 



~l7 



sin(i2f) 



(54) 



which is the well known Mathieu equation |86|. The solutions to this equation have 
an exponentially growing amplitude when Q = 2a)ni- 

There are simple situations in which there is intermode coupling. For instance 
for a cubic cavity of size L, the fundamental mode (1,1,1) is coupled to the mode 
(5,1,1) when the external frequency is 12 = 2(i)\i\. In this case the number of pho- 
tons in each mode grows with a lower rate than that of the uncoupled case OTI 



(55) 



We will describe some additional examples of intermode coupling in the next sub- 
section, in the context of a full electromagnetic model. 

It is worth stressing that the creation of scalar particles in 3D cavities has been 
studied numerically in Ref.|90|. At long times, the numerical results coincide with 
the analytical predictions derived from Eqs.([5T]i and ( [52] i, both in the presence and 
absence of intermode coupling. 



4.2.3 The electromagnetic case 

The previous results have been generalized to the case of the electromagnetic field 
inside a cylindrical cavity, with an arbitrary transversal section f33l. Let us assume 
that the axis of the cavity is along the z-direction, and that the caps are located at 
z = and z = L^{t). All the surfaces are perfect conductors. 
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When studying the electromagnetic field inside these cavities it is convenient to 
express the physical degrees of freedom in terms of the vector potentials A^^^' and 
introduced in Section 2. These vectors can be written in terms of the so 
called "scalar Hertz potentials" as A'™' —zx V(j)^^ and jz/'™' —zx V(j)™. For 
perfect reflectors the b.c. do not mix TE and TM polarizations, and therefore the 
electromagnetic field inside the cavity can be described in terms of these two inde- 
pendent scalar Hertz potentials: no crossed terms appear in Maxwell's Lagrangian 
or Hamiltonian. 

The scalar Hertz potentials satisfy the Klein-Gordon equation. The b.c. of both 
potentials on the static walls of the cavity are 

0™|,=O = O; ^|t,ans -0, (56) 

an 



dij)™ 



,=0 =0; (/)™|trans = 0, (57) 



where d /dn denotes the normal derivative on the transverse boundaries. On the 
other hand, the b.c. on the moving rrdiTor has been already discussed in Section 2 



(see Eqs ( 16 1 and ( 17 1). In terms of the Hertz potentials they read as 



^^\=L-M = ; {d, +4 = 0. (58) 



The energy of the electromagnetic field 



1 



H= — I d\(E^+B^-) = H^^+H™ (59) 



can be written in terms of the scalar potentials as 



87t 



(P) (_ ) (P) + (P)' (- Vi ) (/) (P)' + vi (P) vi (P)' 



(60) 



where dots and primes denote derivatives with respect to time and z respectively. 
The supraindex P corresponds to TE and TM and denotes the gradient on the xy 
plane. 

The quantization procedure has been described in detail in previous papers ll32l 
[33l|9T]. At any given time both scalar Hertz potentials can be expanded in terms of 
an instantaneous basis 

(/.(p)(x,o =i:fl™cr«r(x,f)+c.c., (61) 

n 

where afj^ are bosonic operators that annihilate the IN vacuum state for f < 0, and 

(P) 

C„ are normalization constants that must be appropriately included to obtain the 



usual form of the electromagnetic Hamiltonian (60i in terms of annihilation and 
creation operators. 
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For TE modes, the mode functions are similar to those of the scalar field satisfy- 
ing Dirichlet b.c. described in the previous section 

= LfiL"?TE(0 v/27Z^sin (^g^) v^Jxx). (62) 

For TM modes, the choice of the instantaneous basis is less trivial and has been 
derived in detail in Ref. ||32l 

«r = L[eL"?rM(0+eiTM(0^(^,0]y27z;(0cos r^.ix^). m 

Here the index m ^ is a vector of non-negative integers. The function g{z,t) = 
L{t)L,{t)^{z/Ly{t)) (where ^(z) is a solution to the conditions ^(0) = ^(1) = 
(9,(^(0) = 0, and d,£,{l) — —1) appears when expanding the TM modes in an in- 
stantaneous basis and taking the small e limit. There are many solutions for £,{z), 
but all of them can be shown to lead to the same physical results |32|. The mode 
functions (xx) and (x±), are described below for different types of cavities. 

The mode functions jg^-pM satisfy second order, mode-coupled linear differ- 
ential equations similar to Eq. ( |49] ) ||32ll . As before, for the "parametric resonant 
case" [Q — 2(0j, for some n) there is parametric amplification. Moreover, for some 
particular geometries and sizes of the cavities, different modes n and m can be cou- 
pled, provided either of the resonant coupling conditions Q = \cOn± a)m| are met. 
When intermode coupling occurs it affects the rate of photon creation, typically 
resulting in a reduction of that rate. 

The number of motion-induced photons with a given wavevector n and polariza- 
tion TE or TM can be calculated in terms of the Bogoliubov coefficients. When the 
resonant coupling conditions are not met, the different modes will not be coupled 
during the dynamics. As in the scalar case, the system can be described by a Mathieu 
equation (|54]i for a single mode. As a consequence, the number of motion-induced 
photons in that given mode will grow exponentially. The growth rate is different for 
TE and TM modes 1321 

(-^,te(0> =sinh2(A„TEef) ; {^n^juit)} = sinh^ {X„jMet) , (64) 

where Xn.jE = n^/2(0t, and An.xM = (2<Bn — n?) /2(0„. When both polarizations are 
present, the rate of growth for TM photons is larger than for TE photons, i.e., 
An,TM > An.TE- As in the case of the scalar field, these equations are valid for 
e^hf < 1. 

We describe some specific examples: 

Rectangular section. For a waveguide of length L^{t) and transversal rectangular 
shape (lengths Lx,Ly), the TE mode function is 

v„„,,(x,) = -^cos — j COS — j , (65) 
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with rix and riy non-negative integers that cannot be simultaneously zero. The spec- 
trum is 

(On„ny,n, = \J {n.T: / L,)^ + [uyTZ / Ly)^ + {n,n / L,f , (66) 
with n, > 1 . The TM mode function is 

'"m^,mv(x±) = rj—r- Sin — — sin — — , (67) 



y/ L^Ly \ Lx J \ Ly 

where m^, are positive integers. The spectrum is given by (Om^,my,m,, with m, > 0. 

Let us analyze the particular case of a cubic cavity of size L under the paramet- 
ric resonant condition Q — 2cOk. The fundamental TE mode is doubly degenerate 
((1,0, 1) and (0, 1,1)) and uncoupled to other modes. The number of photons in 
these TE modes grows as exp{7zet^/2L). The fundamental TM mode (1,1,0) has 
the same energy as the fundamental TE mode, and it is coupled to the TM mode 
(1, 1,4). Motion-induced TM photons are produced exponentially as exp(4.4ef/L), 
much faster than TE photons. 

Circular section. For a waveguide with a transversal circular shape of radius R, 
the TE mode function is 

v™ (xx ) - ^ -^===7„ (y„,n § ) (68) 



nm 



where /„ denotes the Bessel function of nth order, and y„m is the mth positive root 
of the equation 7,' {y) = 0. The eigenfrequencies are given by 



-n....-\l[:-f] - (69) 

where n,>l. The TM mode function is 

r™(xx) = ^— ./„(x„,„§) (70) 

y/n RJ„+i{x„,„) \ RJ 



where x„„, is the mth root of the equation J„ (x) — 0. The spectrum is given by Eq.( 69 1 
with y„m replaced by x,„„ and > 0. Denoting the modes by (n, m, n,), the lowest TE 
mode is (1,1,1) and has a frequency COm = [l Ml / R) y^T+2m2{R/L^ . This 
mode is uncoupled to any other modes, and the number of photons in this mode 
grows exponentially in time as exp(7i;ef/^l +0. 343 {L^/R)^Ly) when parametri- 
cally excited. The lowest TM mode (0, 1,0) is also uncoupled and has a frequency 
(Soio = 2AQ5/R. The parametric growth is exp (4.81ef/^). For large enough 
(L- > 2.03/?), the resonance frequency tOm of the lowest TE mode is smaller than 
that for the lowest TM mode. Then the (1,1,1) TE mode is the fundamental oscil- 
lation of the cavity. 
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4.3 Time dependent electromagnetic properties 

From a theoretical point of view, it is possible to create photons from the vacuum 
not only for a cavity with a moving mirror, but also when the electromagnetic prop- 
erties of the walls and/or the media inside the cavity change with time. Given the 
difficulties in a possible experimental verification of the DCE for moving mirrors, 
the consideration of time dependent properties is not only of academic interest, but 
it is also relevant for the analysis of the experimental proposals discussed in Section 
5. 

A setup that has attracted both theoretical and experimental attention is the pos- 
sibility of using short laser pulses in order to produce periodic variations of the 
conductivity of a semiconductor layer placed inside a microwave cavity. The fast 
changes in the conductivity induce a periodic variation in the effective length of the 
cavity, and therefore the creation of photon pairs Il42ll43l . This setup has been an- 
alyzed at the theoretical level [44j|92J|93Jl94|, and there is an ongoing experiment 
aimed at the detection of the motion induced radiation (see Section 5). 

For the sake of clarity we discuss in detail the model of a massless scalar field 
within a rectangular cavity with perfect conducting walls with dimensions L^, Ly, 
and Lj. described in Ref.t44|. At the midpoint of the cavity (x = Lx/2) there is a 
plasma sheet. We model the conductivity properties of such material by a delta- 
potential with a time dependent strength V{t). This is a time dependent general- 
ization of the model introduced in Ref.|66|. The strength of the potential is given 
by 

V(0-4;r^, (71) 

m* 

where e is the electron charge, m* the electron's effective mass in the conduction 
band and n{t) the surface density of carriers. We assume that the irradiation of the 
plasma sheet produces changes in this quantity. The ideal limit of perfect conduc- 
tivity corresponds to V ^ °°, and V — > to a 'transparent' material. The strength 
of the potential varies between a minimum value, Vq, and a maximum Vmax- The 
Lagrangian of the scalar field within the cavity is given by 

^^^-d^^d^<j>-^^5{x-L,/2),p\ (72) 

where 5 (x) is the one-dimensional Dirac delta function. The use of an infinitely thin 
film is justified as long as the width of the slab is much smaller than the wavelengths 
of the relevant electromagnetic modes in the cavity. The corresponding Lagrange 
equation reads, 

{V^-d^)(^=V{t)5ix-L,/2)^. (73) 

We divide the cavity into two regions: region I (0 < jc < Lj/2) and region II 
(Lx/2 < X < Lv). Perfect conductivity at the edges of the cavity imposes Dirichlet 
b.c. for the field. The presence of the plasma sheet introduces a discontinuity in the 
.jc-spatial derivative, while the field itself remains continuous. 
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^\{x = Lx/2,t) = (l)u{x = Lj2,t), 
We will consider a set of solutions that satisfies automatically all b.c. 



V^m(x,f) 



■sin {km,{t)x 




sin 



(74) 



(75) 



where m,. m, are positive integers. The function i/ni depends on t through k„^{t), 
which is the m;f-th positive solution to the following transcendental equation 



2k,„^. tan 



-V{t). 



(76) 



To simplify the notation, in what follows we will write instead of k,n^. Note that, 
when y(f) — > o°, the solutions to this equation become the usual ones for perfect 
reflectors, A;„, — mLx/2, with m a positive integer. 
Let us define 



fm(x,f) 



(77) 



These functions satisfy the b.c. and the orthogonality relations 

('f^m,%) = [1 -sin(^,„(f)L;,)/^,„(f)Lj ^,„. 

There is a second set of modes with a node on the cavity midpoint. As these solutions 
do not "see" the slab, they will be irrelevant in what follows. 

For t <0 the slab is not irradiated, consequently V is independent of time and has 
the value Vq. The modes of the quantum scalar field that satisfy the Klein Gordon 
equation ( 73 1 are 



Mm(x,f) = -^===fm(x,0) , 



(78) 



where (B^ = (k'^)2 + (^^^ + and k^ is the m-th solution to Eq.(76 



for 



y = Vq. At f = the potential starts to change in time and the set of numbers {A:,,,} 
acquires a time dependence through Eq.(76 1. 



Using Eq. ( 78 1 we expand the field operator as 

0(x,O = [^mMra(x,r) +^m«m(x,f)] : 



(79) 



where b„i are annihilation operators. Notice that in the above equation we omitted 
the modes with a node at x = Lx/2 because their dynamics is not affected by the 
presence of the slab. 

For f > we write the expansion of the field mode Ms as 
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«s(x,f >0)=£pi^'(f)fm(x,f). (80) 

m 

Assume a time dependent conductivity given by 

V(f) = ^0 + (Vmax-Vo)/(r) , (81) 

where /(f) is a periodic and non-negative function, f{t) = f{t + T) > 0, that van- 
ishes at f = and attains its maximum at /(Tg) = 1. In each period, f{t) describes 
the excitation and relaxation of the plasma sheet produced by the laser pulse. Typi- 
cally, the characteristic time of excitation Tg is the smallest time scale and satisfies 
^ r. Under certain constraints, large changes in V induce only small variations 



in k through the transcendental relation between k and V (see Eq. (76 1). In this case, 
a perturbative treatment is valid and a linearization of such a relation is appropriate. 
Accordingly we write 

fe„(f)=fe^(l +£„/(?)), (82) 

where 

e„ = "^J^}^^ ^ . (83) 

L,(^o)2 + yo(l + ^) 

The restriction for the validity of the perturbative treatment is VqLx 3> Vmax/Vb > 1- 
These conditions are satisfied for realistic values of L^, Vq, and Vmax- 

Replacing Eq.(80i into (V^ — d})us = we find a set of coupled differential 

is) 

equations for the amplitudes Pm (f). The dynamics is described by a set of coupled 
harmonic oscillators with periodic frequencies and couplings, as already discussed 
in this section. It is of the same form as the equations that describe the modes of a 
scalar field in a three dimensional cavity with an oscillating boundary. For the same 
reasons as before, a naive perturbative solution of previous equations in powers of £„ 
breaks down after a short amount of time when the external frequency is tuned with 
some of the eigenfrequencies of the cavity. Assuming that /(f) is a sum of harmonic 
functions of frequencies Qj = j2n/T, the resonance condition is Qj = 2S)a for 
some i and n. If there is no intermode coupling, a nonperturbative solution gives an 
exponential number of created photons in that particular mode 

(^n(O) = {blbn) = '£2a)n\A^\t)\^ « sinh2 (^^|^e«f) , (84) 

where /,■ is the amplitude of the oscillations of /(f) with frequency Qj. 

A full electromagnetic calculation has been presented in Ref. ||94l . It was shown 
there that the scalar model presented here describes the TE electromagnetic modes 
inside the cavity. The treatment of TM modes involves an independent scalar field, 
with a potential proportional to 5'(x — L,/2). Moreover, the model has also been 
generalized to the case of arbitrary positions of the plasma sheet within the cavity 
The number of created TE photons depends strongly on the position of the 
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layer, and the maximum number is attained when it is located at the midpoint of the 
cavity. On the other hand, for TM modes this dependence is rather weak. 

In the treatment above no dissipation effects are considered (the delta-potential 
is real). Similar calculations for lossless dielectric slabs with time-dependent and 
real permittivities Il93l also neglect dissipation. However, it has been pointed out 
that dissipative effects may be relevant in the evaluation of created photons 1921 . 
In general, one expects the electromagnetic energy to be dissipated in the cavity 
walls, in the plasma sheet, and/or in dielectric slabs contained in the cavity. In res- 
onant situations without dissipation, we have seen that the dynamics of the relevant 
electromagnetic mode is described by a harmonic oscillator with time dependent 
frequency. A phenomenological way of taking into account dissipative effects is to 
replace this equation with that of a damped oscillator |95| . Of course this model 
cannot be consistently quantized unless one includes a noise term, otherwise the 
usual commutation relations are violated. Using the quantum noise operator ap- 
proach |96 1 one can estimate the rate of photon creation in this model and, provided 
the dissipation is not too large, the number of photons still grows exponentially, 
although at a smaller rate. However, it has been recently argued |38| that these re- 
sults should be valid only in the short time limit, while in the long time limit the 
system should reach a stationary state with a constant number of photons inside the 
cavity. As the calculations in li38l involve ID cavities, this point deserves further 
investigation. 



5 Experimental perspectives 

Since the first theoretical predictions about motion induced radiation, it was clear 
that the experimental observation of this effect was not an easy task. As mentioned 
at the end of Section 2, the photon creation produced by a single accelerated mirror 



is extremely small in realistic situations (see Eq. (20l). 

The most promising situation seems to be the photon creation by parametric am- 
plification described in Section 4. However simple numerical estimations show that, 
even in the most favorable cases, it is difficult to observe the DCE in the laboratory. 
In all the 3D examples discussed in Section 4, the number of created photons grows 
exponentially in time as 

(^) = sinh2(Tja)ef) , (85) 

where O) is the frequency of the resonant mode and T] is a number of order 1 re- 
lated to the geometry of the cavity. Here £ denotes the relative amplitude of the 
oscillations in the moving mirror case, or the relative amplitude of the oscillations 
of the relevant component of the wavevector in the case of time dependent conduc- 



tivity (see Eq.(|82]i). This equation is valid as long as e^COt <C 1 and neglects any 
dissipative effects. As the electromagnetic cavity has a finite g-f^ctor, a rough es- 
timation of the maximum number of created photons {.yVnwx) is obtained by setting 
tmax = 2/ft) in the above equation. As mentioned at the end of Section 4, there is 
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no agreement in the literature about this estimation. Calculations based on the use 
of a master equation ||39l give an exponential growth with a rate diminished by a 
factor r = 1 — 1/ {2Qe) (see also Ref. |'97l). On the other hand, it was shown that, 
in the case of ID cavities, the total number of photons inside the cavity should reach 
a constant value proportional to the finesse of the cavity at long times 1381 . It was 
argued in the same work that the exponential growth in the presence of dissipation 
would be valid only at short times. In any case, it is clear that a necessary condition 
to have an observable number of photons is that 2Qe > 1. 

Assuming a cavity of length Lq ~ 1 cm, the oscillation frequency of the mirror 
should be in the GHz range in order to meet the parametric resonance condition. 
A plausible possibility for reaching such high mechanical oscillation frequencies is 
to consider surface vibrations, instead of a global motion of the mirror |30|. In this 
context, the maximum attainable values of the relative amplitude would be around 
e ~ 10^^, and therefore the quality factor of the cavity should be greater than 10^ 
in order to have a non-negligible number of photons. Microwave superconducting 
cavities with Q-factors as high as lO'^ have been built |98|. However, the Q-factor 
would be severely limited by the presence of an oscillating wall. Therefore, it is 
an extraordinary challenge to produce extremely fast oscillations while keeping the 
extremely high Q-factors needed in the DCE. Moreover, the oscillations should be 
tuned with high precision to parametric resonance with a cavity mode. 



5.1 High frequency resonators and photon detection via 
superradiance 

A concrete setup for producing and detecting motion induced photons has been 
proposed in Ref |41|. A Film Bulk Acoustic Resonator (FBAR) is a device that 
consists of a piezoelectric film sandwiched between two electrodes. An aluminum 
nitride FBAR of thickness corresponding to a half of the acoustic wavelength can 
be made to vibrate up to a frequency of 3GHz, with an amplitude of £ = 10^**. The 
expected maximum power of Casimir photons produced by such a FBAR depends 
of course on the Q-factor of the cavity. It can be estimated to be 

Assuming that Qe = 0(1), this gives P„,„, ~ lO^^^W, which is too small for direct 
detection. 

However, this low power could be detected using ultracold atoms. Let us consider 
a cavity filled with an ensemble of population-inverted atoms in a hyperfine state 
whose transition frequency is equal to the resonance frequency of the cavity. Then 
the Casimir photons can trigger a stimulated emission of the atoms, and therefore 
they can be indirectly detected by this form of superradiance. Ref. BTI contains a 
description of the experimental setup that could be used to observe the DCE, and 
a discussion about the rejection of signals not produced by the Casimir photons. 
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In particular, stimulated amplification could also be triggered by the spontaneous 
decay of one of the atoms (superfluorescence). In order to discriminate between 
both effects it could be necessary to attain larger values of Qe. 



5.2 Time dependent conductivity induced by ultra-short laser 
pulses 

In order to avoid the experimental complications associated with the high frequency 
motion of the mirror, it is possible to produce effective changes in the length of 
the cavity by inducing abrupt changes in the reflectivity of a slab contained in the 
cavity, as already mentioned in Section 4.3. This can be done by illuminating a 
semiconducting slab with ultra short laser pulses B2ll43]| . An experiment based on 
this idea is currently being carried out by the group of Padova 1*991. 

In this case, a numerical estimation of the maximum number of photons created 
in the cavity looks, at first sight, much more promising than in the case of moving 
mirrors. Using a slab of a thickness around 1 mm, it is possible to reach values of 
£ as large as 10^^, and therefore the constraints on the g-f^ctor of the cavity are 
considerably milder. Moreover, it is experimentally possible to generate trains of 
thousands of laser pulses with a repetition frequency on the order of a GHz. 

In Padova's setup (see Fig. 4), a high Q « 10^ superconducting cavity contains 
a GaAs semiconducting slab. The laser pulses are tuned at 4.70 GHz, twice the fre- 
quency of the fundamental TE mode of the cavity. However, it has been pointed 
out ||92ll that dissipative effects may play an important role in this kind of experi- 
ment. Indeed, the changes in the conductivity of the slab are due to the creation of 
electron-hole pairs by the laser pulses, and during this process the dielectric permit- 
tivity acquires an imaginary part. The associated dissipation prevents photon cre- 
ation unless severe constraints on the properties of the semiconductor are fulfilled: 
it must have a very short recombination time (tenths of ps), and a high mobility 
(around lm^{Vs)^^). A slab with these characteristics has been constructed by irra- 
diating a GaAs sample with fast neutrons, in order to reduce the recombination time 
of the original sample (about 1 ns) while keeping constant the value of the mobil- 
ity [99] . Photons are detected using a loop antenna inside the cavity. The minimum 
number of photons that can be detected is around 100, below the expected signal of 
Casimir photons |l47^ '99l. 

A related setup is illumination of a superconductor instead of a semiconductor 
surface. The advantage in this case is that dissipative effects are less important, be- 
cause the variation of the imaginary part of the permittivity is much smaller for 
superconductors than for semiconductors in the microwave region 1 100 1. Moreover, 
since the abrupt changes in the conductivity are due to local heating of the surface 
(and not to the creation of electron-hole pairs as in the semiconductor), the inten- 
sity of the laser can be considerably smaller, reducing unwanted effects of energy 
accumulation inside the cavity. 



32 



Diego A. R. Dalvit, Paulo A. Maia Neto, and Francisco Diego Mazzitelli 



5.3 Optical parametric oscillators 

Standard nonlinear optics can be interpreted, in some cases, as a time-dependent 
modulation of the refractive index. Ref. |38| considered an optical parametric os- 
cillator (OPO) with a pump laser beam of frequency Q and amplitude fipump in- 
teracting with a very thin x^^^ nonlinear crystal slab placed on the interior side of 
the cavity mirror. For a type-I arrangement, the total polarization component along 
a suitable crystal symmetry axis may be written in terms of the intracavity electric 
field component along the same direction as ll38l 

^(0 = eo (^Z*'' + ^z'^*£pumpsin(i2f - 9)^ E{t), (87) 

where eo is the vacuum permittivity, j''' and x^^^ are the relevant components of 
the linear and second-order nonlinear susceptibility tensors, and is a phase that 
depends on the pump beam phase and the position of the crystal. 




Fig. 4 Superconducting cavity of the experimental setup of MIR experiment at Padova to measure 
the dynamical Casimir effect (Courtesy of Giuseppe Ruoso). 
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The total susceptibility as given by eq. ( [ST] ) (including the nonlinear second-order 
term) corresponds to an effective refractive index oscillating at the pump beam fre- 
quency 12 , thus leading to a modulation of the optical cavity length. This is formally 
equivalent to modulating the physical cavity length by bouncing the mirror with fre- 
quency Q.. But in the OPO the pump beam frequency is in the optical range, as 
are the generated photons with frequencies satisfying o + o' ^ Q. For (O and (o' 
corresponding to cavity modes, parametric amplification is enhanced and the result- 
ing photon flux is typically several orders of magnitude higher than in the case of 
mechanical motion [38 1. 



5.4 Superconducting coplanar waveguides 

Another possibility to induce fast variations of the b.c. on the electromagnetic field 
is to consider a coplanar waveguide terminated by a SQUID B9l |5P|. A time- 
dependent magnetic flux can be applied to control the effective inductance of the 
SQUID, which in turn produces a time-dependent Robin boundary condition for the 
phase field (time integral of the electric field), equivalent to that of a transmission 
line with a variable length. This setup simulates a moving Robin mirror in ID with 
an effective velocity that might be close to the speed of light. As a consequence, 
the first-order non-relativistic results | l22l l64l . based on the perturbative approach 
outlined in Section 2, must be modified by the inclusion of higher-order frequency 
sidebands |49|. 

As in the previous examples, it is crucial to check if the flux of Casimir photons 
can be discriminated from other sources of photons, like the classical thermal con- 
tribution. The analysis presented in Ref.|49| shows that this is the case, for realistic 
values of the parameters, at temperatures below 70 mK. 



6 Final remarks 

We have reviewed some theoretical and experimental advances in the analysis of 
moving bodies or time dependent boundary conditions coupled to the vacuum fluc- 
tuations of the electromagnetic field. 

Accelerated neutral bodies produce the emission of real photons, while experi- 
encing a radiation reaction force. When the dynamics of the bodies is treated quan- 
tum mechanically, the interaction with the vacuum fluctuations not only causes this 
dissipative force, but also an appreciable amount of decoherence, which is a con- 
sequence of the entanglement between the mirrors and the electromagnetic field. 
This is a particular example of quantum Brownian motion, where the Brownian par- 
ticle (mirror) loses coherence while being subjected to a damping force due to its 
coupUng to the environment (the quantum field). 
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When two neutral bodies are in relative motion, we expect velocity-dependent 
forces between them. There is a particularly interesting situation in which two par- 
allel, non-perfectly conducting slabs are in relative parallel motion with constant 
velocity. In this case, there is a vacuum friction between the slabs even in the ab- 
sence of real photons. The effect can be understood in terms of the interaction of 
image charges, or as the interchange of virtual photons between the surfaces. There 
are similar friction forces for neutral atoms moving near surfaces. 

The rate of photon creation produced by a single accelerated body in free space 
is deceptively small in reahstic situations. However, in closed cavities a much larger 
number of photons may be produced by parametric amplification. Indeed in an ideal 
3D cavity one expects an exponential growth in the number of photons when its size 
varies periodically at an appropriate resonant frequency, making detection of photon 
creation not an impossible task. The calculation of photon creation in the presence 
of ideal conductors have been performed in ID and 3D using different analytical 
approximations. The results are consistent and have been confirmed by fully niuner- 
ical calculations. However, the case of moving mirrors with finite conductivity (i.e 
electromagnetic cavities with a finite Q-factor) is not a completely settled issue. In 
ID cavities, at short times the growth of the total number of created photons is still 
exponential (with a different rate), while at large times the total number of photons 
should reach saturation. This problem has not yet been solved for 3D cavities. The 
difficulties in evaluating the DCE for mirrors with finite conductivity recalls a sim- 
ilar situation in the static Casimir effect, where the evaluation of the Casimir force 
depends strongly on the theoretical model used to describe the conductivity of the 
bodies, and there are interesting correlations between finite conductivity, tempera- 
ture and geometry. These correlations may have relevant counterparts in the dynam- 
ical problem. In any case, although difficult, the direct experimental detection of the 
motion induced radiation is not out of reach, as long as one can keep a very high 
2-factor in a cavity with moving walls. In particular, there is a specific proposal that 
involves nanoresonators in a high Q-cavity filled with a gas of cold atoms to detect 
a small number of photons through superradiance. 

An exponentially large number of photons can also be produced when some elec- 
tromagnetic property of the cavity varies periodically with time. Of particular im- 
portance is the case in which the conductivity of a semiconductor or superconductor 
slab placed inside an electromagnetic cavity is modulated using short laser pulses. 
Theoretical estimates show that this setup could be implemented, with milder re- 
quirements on the 2-factor of the cavity. Once again, there is no comprehensive 
theoretical model that takes into account the (dissipative) response of the slab to the 
laser pulses, and its relevance for the photon creation process. However, this is a 
promising alternative and there is an ongoing experiment at Padova based on this 
setup. 

There are other possibilities to produce fast variations of the boundary condi- 
tions on the electromagnetic field, that involve optical parametric oscillators or su- 
perconducting waveguides. The theoretical analyses suggest that it should be easier 
to detect the photons created in these setups than in the case of moving mirrors. 
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In summary, there is a plethora of interesting effects related to the electromag- 
netic vacuum fluctuations in the presence of moving bodies and/or other time depen- 
dent external conditions. The eventual experimental confirmation of some of these 
effects will certainly produce an increasing activity on this subject in the near future, 
as was the case for the static Casimir effect following the first realization of precise 
experiments in that area since 1997. 
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